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I. INTRODUCTION. 



In the two decades following its invention, scanning tunnelling microscopy (STM) has 
proved to be a valuable tool for investigating surfaces on an atomic scale. More recently, 
several experiments show a growing interest in the study of structures that are situated 
in the bulk below the surface in both semiconductors and metals. Whereas in the former 
(i.e. semiconductor) case the absence of effective screening allows dopants down to the third 
subsurface layer to be viewed directly as apparent topographic features, [1] the situation for 
metals turns out to be somewhat more complicated. One method that has been suggested for 
imaging structures buried in metal involves several surface study techniques to be employed 
simultaneously in combination with STM [2]. However, although this experiment has lead to 
successful identification of subsurface defects, it cannot be used as a tool for probing the exact 
depth. Also, its employabihty is hmited to certain specific alloys only. A more successful 
approach however seems to be by probing standing electron waves[3, 4]. The groundwork 
of these experiments is as described in Ref. [5], where Cu(lll) surface states form a two- 
dimensional nearly free electron gas. When scattered from step edges or adatoms, these 
states then form standing waves which can be probed by scanning tunnelling spectroscopy 
(STS). 

Although it has been proposed to utilize these surface states for imaging subsurface 
impurities [6], the exponential decay of the wave function amplitudes into the bulk will limit 
the effective range to the topmost layers only. Bulk states however, of which the square 
falls of with only r^, form a good alternative. To demonstrate this, we mention results that 
were obtained by bulk state spectroscopy on relatively large structures such as Ar bubbles 
submerged in Al. [3] and Si(lll) step edges buried under a thin film of Pb [4]. In these 
experiments, bulk electrons are found to be confined in a vertical quantum well between the 
surface and the top plane of the object of interest. 

In this paper we show that the investigation of the nonlinear conductance of a point 
contact placed on a metal surface makes it possible to determine the position of point- 
like defects such as the vacancies or foreign atoms inside the metal in the vicinity of the 
contact. We consider theoretical models both for the cases of a tunnel point contact and for 
a baUistic quantum contact. We look for conductance oscillations caused by interference of 
electrons that arc transmitted directly, and electrons that are first backscattered elastically 
by the defect and again scattered forward by the contact (i.e. the tip-sample junction), 
much in the same way as was described for atomic point contacts in Refs. [7-9] The effect of 
such quantum interference on the nonlinear conductance of quantum wires was theoretically 
analyzed in Refs. [11-13], but the point contact geometry was not studied yet. 

The organization of this paper is as follows. In Sec. II we consider a tunnel junction 
in the limit of a high potential barrier. The interaction of the transmitted electrons with 
a single impurity near the junction is taken into account by perturbation theory with the 
electron-impurity interaction as the small parameter. A general analytical expression for 
the voltage dependence of conductance, G iV). is obtained. It defines G (V) in terms of the 
contact diameter, the distance between contact and the impurity and the parameters that 
characterize the metal, and the transmission of the tunnel junction. In Sec. Ill the conduc- 
tance of a ballistic quantum contact of adiabatic shape is analyzed. In absence of a barrier 
inside the contact electrons can still be reflected from it due to the variation of the confining 
potential. The infiuence on electron scattering by an explicit barrier potential in the center 
of the contact is also discussed. As in Sec. II assuming the electron-impurity interaction 
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to be small wc derive an expression for G iV) and its dependence on the position of the 
defect. In Sec. IV we conclude by discussing the possibilities for experimental exploitation 
of the conductance fluctuations for sub-surface imaging as well as the technical difficulties 
involved. 



II. TUNNEL POINT-CONTACT. 



Let us consider as a first model of our system a nontransparent interface located at 2; = 
between two metal half-spaces, in which there is an orifice (contact), as illustrated in Fig. 1. 
The potential barrier in the plane 2; = is taken to be a 5— function: 

U (r) = Uf (p) 5 (z) , (1) 

where p = {x,y) is a two dimensional vector. The function / (p) — cxd in all points of the 
plane z — except in the contact, where / (p) = 1. At a point r = ro in vicinity of the 
interface, in the half-space 2; > 0, a point-hke defect is placed, see Fig. 1. The electron 
interaction with the defect is described by the potential 

gir)=gS{r-ro), (2) 

where g is the constant of the electron-impurity interaction. In this section we consider 
tunnel junctions and assume that the transmission probability of electrons through the 
orifice is small. In that case the applied voltage drops entirely over the barrier and we 
choose the electric potential as a step function V (z) — V Q {—z) , and take eV > 0. 

The electrical current / (V) can be evaluated [10] from the electron wave functions of the 
system, 

/ {V) = JdkJdS Im (^:V^k) {Q [riF {e k) - np (^k + eV)] . (3) 

Here, Sk = ^^k^/2m* is the electron energy, k is the electron wave vector and m* is an 
effective mass of electron; np (^k) is the Fermi distribution function. The real space inte- 
gration is performed over a surface overlapping the contacts in the region 2; > 0. At low 
temperatures the tunnel current is due to those electrons in the half-space 2; < having an 
energy between the Fermi energy, Sp, and Sp + eV, because on the other side of the barrier 
only states with Sk ^ are available. 

The wave function ip^^ satisfies the Schrodinger equation 

2m* 

v'V'k (P, z) + ^[e^-U (r) - ^ (r) - eV (z)] V'k (p, z) = (4) 

where the wave vector k = {x,kz) has components x and kz parallel and perpendicular to 
interface, respectively. As shown in Ref.[14], Eq. (4) can be solved for arbitrary form of the 
function / (p) in the limit 1/U — > 0. The wave function 'il)-^^{p,z) for A;^ > in the main 
approximation of the small parameter ~ 1/C/ takes the form: 

^i, (p, z) = e-" {f^ - e-'~'^) + i^-) (p, z) {z<0), (5) 

^|^y.{p,z) ^ ^cpi^^p^z) {z>0\ (6) 
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FIG. 1: Model of a tunnel junetion eontact as an orifice in an interface that is nontransparent for 
electrons except for a circular hole, where tunnelling is allowed. Trajectories are shown schemat- 
ically for electrons that are reflected from, or transmitted through the contact and then reflected 
from a defect. 



where k= ^x, —k^ , k — y/k^ + 2meV/h^. The function ip^^ (p, z) satisfies the conditions of 
continuity and the condition of the jump of its derivative at the boundary z — 0. At large 
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U these conditions are reduced to 



(/^' 0) = ^k^^ (P, 0) ; (7) 



m 



ik = ^/ (p) V^i^^ (p, 0) . (8) 
In the absence of the defect {g — 0) the wave function was obtained in Ref.[14], 

^ oo 
— oo 

where 

oo 

/ i(>f-V)p 
rfp-^^, (10) 

— oo 

and A;^ = a/ + A;^ — x'^ . For a circular contact of a radius a, defined by / (|p| ^ a) = 1 
and / (IpI > a) — >■ oo, the function F — x') takes the form 

^^^_^^^2aJ^ax-x>)_ (11) 
|x — x'l 

In order to introduce the effect of the impurity we solve the Schrodinger equation for the 
Fourier components $k (x^, z) of the function ip^^ {p, z) 

oo 

<^L^^(P,^) = / cixe-^^k (x, ^) . (12) 

— oo 

For z > this equation takes the form 

-x'2$k (x', 2;) + + [^k*k ( x', ^) - (13) 



^5(z-^o)e-^^^vL"'^(Po,^o) 



= 



Integrating Eq. (13) near the point z = ZQwe obtain the effective boundary condition: 

(x', z ^ + 0) - (x', z = - 0) = ^e-^^o^^+) (p^, (14) 

To proceed with further calculations we assume that the electron-impurity interaction con- 
stant g is small and use perturbation theory. In this approximation we replace ip^^'^ by 
'^ok'' (9) ill the right hand side of Eq. (14). Solving the Schrodinger equation (13) with the 
boundary conditions (7), (8), (14), and the condition of continuity of the function $k (x, z) 
&X z — Zq, we obtain in the region z > Zq 

$k(x',^) =^k(x')e^'=^^ (15) 

where _ 

^ (X - x') - ^^(^^ (po, .0) e-^-'-o sin (A;;.o) ■ (16) 
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Using Eq. (13) we find the wave function 



oo 



(17) 



where Zi = z, Z2 = Zq, when z > Zq, and Zi = Zq, Z2 = z when < 2; < ^o- The 
modulus of the wave function (17) in a plane through the impurity and normal to the 
circular contact is illustrated in Fig. 2, for an incident wave vector normal to the interface 
{kz = kp = j^^y2m*eF)- One recognizes a interference pattern of partial waves reflected at 
the impurity with those emanating from the contact. These contain the information that 
we hope to extract. 

Substituting wave function (6) into Eq. (3) and taking into account Eq. (17), we calculate 

the current-voltage characteristics I{V). After integration over all directions of the wave 
vector k and integration over the space coordinate p in a plane z = const {z > zq), retaining 
only terms to first order in g (i.e. ignoring multiple scattering at the impurity site), the 
current is given by 



I{V) 



2eh^ 



00 

/ 



dkk^ 



dpidp2 



f (Pi) / (P2) 



[riF (^k) - riF (^k + eV)] 



(18) 



A'^ikp) ^_m*gkA{kXi)A{kX2) ^2 f^ikp) , 2m*ey sin (A;p) 

kp 



-27r 



+ 



where p = jp^ - p^\ , Ai = ^/zf+~\p^^^p^ , X2 = ^^0 + |Po -P2f, and 



A{x) 



smx 



— cosx 



X 



(19) 



Differentiating Eq. (18) with voltage V and integrating over the absolute value k of the wave 
vector, in the limit of low temperatures, T = 0, we obtain the conductance G (V) of the 
system 



GiV) 



e^klh^ 



dp^dp2 



n{m*UfJJ /(Pi)/(p2) 



A' 



(kFp) 



(20) 



27r- 



m*gkF ^ ) ^ (,^^^2 j 2 M r"^^ ) 2m*e1/ y^^P 



X^X^ 



8nm*g 2 



P' 



+ 



^2 



kpp 



k^dk 



sin {kp) A (kXi) A {kX2) 



where kp = \/kp + 2m*eV/h? is the Fermi wave vector accelerated by the potential differ- 
ence, and we have assumed eV/sp < 1. For Aj and p much larger than the Fermi wave length 
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FIG. 2: Modulus of the wave function in the vicinity of a tunnelhng point-contact in a plane 
perpendicular to the contact axis, having an impurity at (pg = 5, 20 = 15). The incident wave has 
a wave vector normal to the point contact a = A_p/47r. 

\f = ^/kp, the function A oscillates with {eV/eF)kF)^i and {eV/eF)kFp, which results in 
a fluctuation of the conductance with applied voltage. The first term is square brackets of 
Eq.(20) describes the conductance G (V) = Gc iV) in the absence of a defect ((7 = 0). 

For a contact of small diameter a <^ \f Eq. (20) can be simplified and the conductance 
is given by 
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FIG. 3: Dependence of the normalized conductance G/Gc for a single tunnel point contact as a 
function of the position of the defect (/Sq, -^o)) contact radius is a = Ai;'/47r . 



G{V) = G,{V){1 



6nm*g 



) 




Qm*eV 



(21) 



12'Kgm* zl 




where 7 = y -Zq + IPoI ^^e distance between the contact and the defect. The conductance 
Gc {V) of the tunnel junction of small cross section S = no? -C is given by 



for small transmission coefficient (k) — ti^k'^ / m*'^U'^ <IC 1. 

In numerical calculations we use a value for the dimensionless parameter 2Tcm* gk p / ti? — 
0.3 to characterize the strength of the defect scattering. Fig. 3 shows a plot of the dependence 
of the normalized conductance G (0) /Gc (0), Eq. (20), for the contact as a function of the 
position of the defect (pg, 2^0) in the hmit of low voltage F — > 0. We observe a suppression 
of the conductance that is largest when the contact is placed directly above the defect and 
find that G is an oscillatory function of the defect position. In Fig. 5 we show the voltage 
dependence of the normalized conductance G (V) /GciV). Eq. (20), for Pq = and as a 
function of the depth zq of the defect under the metal surface. 



GeiV) 



9hU^ 




(22) 
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FIG. 4: 




FIG. 5: Voltage dependence of the normalized conductance G {V) /Gc{V) of a tunnel point contact 
for Pq = and as a function of the depth zq of the defect under metal surface, contact radius is 
a = Xf/Att. 

III. BALLISTIC CONTACT. 

In this section we consider another limit of a junction, a cylindrically symmetric, ballistic 
contact of adiabatic shape, Fig. 6. The center of the contact is characterized by a (5-function 
potential barrier of amplitude U. In one of the banks of the contact a single defect is situated 
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FIG. 6: Model of a ballistic contact of adiabatic shape having a defect sitting nearby. The tra- 
jectories of electrons that are transmitted through the contact and reflected from the defect are 
shown schematically. 

at the point fq = (po- ^o) - in the half-space Zo > 0, such that the distance 7 between the 
center of the contact r = and the defect is much larger than the characteristic length L 
of the constriction (see Fig. 6). The shape of the contact is described by the radius as a 
function of the z-coordinate, a {z) . The contact size is given by a (0) = Gq, while a{z) —>■ 00 
for 1^1 — >■ 00. The adiabatic condition implies that the radius of the contact a{z) varies 
slowly on the scale of the Fermi wavelength. As a result, the electric potential V (r) drops 
dominantly over the same characteristic length L, as can be derived from the condition of 
electroneutrality. In the Landauer formalism the exact distribution of V (r) is not important 
for determining the conductance of a quantum constriction, which can be expressed using 
only the difference of potentials V in the banks far from the contact. We will consider 
the effect of quantum interference on the conductance under conditions eV/sp ^ 1 and 
kFjieV/sp) > 1- Fluctuations of G {V) results from the phase shift Aip that the wave 
function accumulates after being scattered by the defect and reflected by the contact, Aip 
~ {eV/eF)kF'y- If 7 3> L, the main part of the electron trajectory is situated in the region 
where the local electric potential V (r) differs only little from its value V = in the bank at 
z ^ 00, and we neglect this small variation of the potential. Assuming hard wall boundary 
conditions, we need to solve the Schrodinger equation, 

2m* 

V^V'« (P, z) + ^[e-g5iv- Vo) - U5 (z)] (p, z) = 0, (23) 
with the boundary conditions 

^^{\p\=a{z)-z) = Q, (24) 

and OL represents the full set of quantum numbers. 

In the adiabatic approximation [15, 16] the "fast" transverse and "slow" longitudinal 
variables in Eq. (23) can be separated and the wave function takes on the form 

V'a (P, = (p, Z) {z) , (25) 
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where /3 — (m, n) is a set of two discrete quantum numbers, which define the transverse local 
eigenvalues Ef^ {z) and eigenfunctions '(/'_|_^ (p, z) . The function ipj_f^ (p, z) depends on the 
coordinate 2; as a local parameter, and its derivatives with respect to z are small. Therefore 
Eq. (23 ) can be separated into two equations. 



2m* 

VpV'±^ (p) = (a) ; 

c?V/3e ^ 2m* 

The functions -^j^^ (p) and 99^^ (2;) satisfy the following conditions: 



^±/3(|p| = a) = 0; 



(^) 



20+0 



,.mc /„\ 



_ 2m*5f 

20-0 " 

for 2; - 



k±/3 (Po)|V/3£(^o); 



-00; 



C^V'fle (^) 



+0 



dz 



2m* U 



(0) ; 



(26) 
(27) 

(28) 

(29) 
(30) 
(31) 



where k — \/2m*e/h and ^p"^^ [z] is the incident wave. Condition (30) means that we 
consider a wave </7^'^ {z) of unit amplitude, which moves from —00 towards the contact. 

For the subsequent calculations we make the explicit choice for the shape of the contact 
a{z) = aocosh(2;/L) . The condition of adiabaticity for this dependence a{z) is L ^ Xp. 
The solution of Eq. (26) is given by 



{z) Jm+i (7 mnJ 



It. 



P 



a yz) 



having eigenvalues 



2^,2 

mn 



2m*a? [z) ' 



n = 0,l,2,...; m = 0,±l,±2. 



(32) 



(33) 



Here, we use cylindrical coordinates p — (p, z) and is the n—th zero of Bessel function 
Jjn- The energy spectrum (33) describes the quantized energy levels inside the constriction 
{z < L) and a quasi-continuous spectrum at 2; ^ L (the distance between the levels Ae^ 



at \z\ 



oc). 



First, we consider a contact without an explicit barrier ([/ = 0). A general solution for 
the longitudinal wave function ip^^ {z) in Eq. (27) has the form. 



<pp,{z) ^ A{i-ey 
B{i-ey 



ikL 

- 2 p 



-ikL — s; —ikL + s + 1; —ikL + s; 



+ 



(34) 



ikL / l — ^ 
1 ^ 



ikL 



F -s;s + l;l + iA;L; 



where F[a,h,c]^) is the hyper geometric function, ^ = tanh (|^) , and s = 
1^—1 -\-i-\J (2L7^„/ao)^ — 1^ . The constants A and B can be found from the conditions 
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(29) and (30). By using the asymptotic form of hypcrgcomctric function at z > zq ^ L and 
in the hmit of a small electron- impurity interaction constant rrfghp/f^ -C 1, we find 



Akz. 



(35) 



where rp and tp are the amplitudes of the reflected and transmitted waves far from the 
contact 



r {-ikL -s)T {-ikL + s + 1) r (ikL) ^ 

r(-s)r(s + i)r(-ikL) 

r {-ikL - s) r {-ikL + s + 1) 



r {-ikL) r (1 - ikL) 

The expression for the current takes the form: 

e 



(36) 
(37) 

(38) 



Substituting the wave function ip^{p,z) from (25) into Eq. (38), with ip±p{p,z) given by 
(32), and (pg^ (z) by (35), and carrying out the integration over p ^ a{z) ai z > Zo, we find 
the total current trough the contact for kL ^ 1, and L ^ ao, 



HV)-J J dsY^TpW {e)-nF {e + eV)] 



X 



(39) 



1 + 



gm 



\i^^l3 (Po)r V(}\ COS (2fc^o + 



Here, 



Tp{e) 
\rp{e)\ 



1 



1 + exp [27rL {kp - k)] ' 
1 

•>/exp [2-kL {k - kf})] + 1 ' 



(40) 
(41) 



^p {e) = {k + kg) L {1 - In ((fc + kp) L)} (42) 
- (A; - A;^) LV^ - 2A;L(1 - lnA;L); 

kp = ^J2/mFFpJfi) / h is the quantized momentum of the transverse electron motion at the 
centre of the contact; i) {\) = V (i) /T (|) ^ -1.96. 

The conductance G {V) = dl/dV, in the low temperature limit, is given by 

G {V) ^GoJ^Tp (sf) { 1 + (Po) r cos (2^^^o + Vp)Y (43) 

where Gq = 2e^ /h is the quantum of conductance, and kp = \/kp + 2meV/h?. All energy 
dependent functions are taken at e — Sp. We also used the condition eV/ep <S 1- 
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FIG. 7: Dependence of the oscihatory part of normahzed conductance AG (0) jGc (0) of the bahistic 
point contact with the position coordinates of the defect (po) -^o)- 

The transmission coefficient is exponentially small [e) ~ exp [— 27rL ik^ — A;)] for e <^ 
€p (0) , while Tp (e) 1 above this energy, e > (0) . For \e - (0)| < (0) Eq. (40) 
agrees with the formula for the transmission coefficient that can be obtained in such case for 
an arbitrary dependence a {z) , by the using an expansion near the point of minimum cross 
section, a (0) (see, e.g., Ref. [18]). For very long constrictions, L 00, Eq. (40) transforms 
to a step function 6 (e — (0)) . For large L the electrons are strongly reflected by the 
contact when k ~ kp. Hence, for observation of conductance oscillations in an adiabatic 
ballistic constriction the contact diameter should be chosen in such a way that Sp ^ (0) , 
i.e., not very far from the middle of a conductance step. 

In the case ?7 7^ the boundary condition (31) must be taken into account. At k ^ kjs 
reflection due to the shape of the contact is negligibly small, as discuss above, and the 
conductance G (V) is described by the same equation (43), but with 

Tp (e) = ^ o, (44) 

1+ {m*U/h^k'f^) 

nfUlh^k' 

\r,{e)\ = ^— ' ' (45) 



y^l + {m*U/h?k'p) 
^(3 (^) = arcsin ( I , (46) 




where fc^ = ^k^ — k^. 

Figure 7 shows the dependence of the oscihatory part of conductance AG{V) = G {V) — 
Gc iy)^ Eq. (43), on the position of the defect at low voltage, — 0, for a contact without 
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ag/g 




FIG. 8: Voltage dependence of the oscillatory part of normalized conductance AG {V) /Gc {V) for 
an adiabatic ballistic point contact for po = as a function of the depth zq of the defect under 
metal surface. 



barrier {U = 0). Here, Gc = Gq ^pTp is the conductance in the absence of a defect {g = Q). 
Figure 8 shows the dependence of the AG(V) on applied bias voltage for a defect sitting on 
the axis of the contact (pg = 0) , and as a function of the distance Zq from the contact center. 
In creating the plots of Figs. 7 and 8 we have used dimensionless parameter 277111* gkp/h'^ = 
0.5 , 27rao = 2.405Ai7' and 2-kL — lOXp, corresponding to a contact having one allowed 
quantum conductance mode. 



IV. DISCUSSION. 



The presence of an elastic scattering center located inside the bulk, either in the vicinity 
of a tunnel contact in an STM configuration or in one of the banks of a ballistic point contact, 
has been shown to cause oscillatory fluctuations in the conductance of the junction. For small 
contact radii (a <^ Xp), these oscillations result solely from interference of electron waves 
that are directly transmitted on the one hand, and electrons that are both backscattered 
by the defect and again reflected by the contact on the other. What now follows is a 
discussion whether this effect can be employed experimentally for three dimensional mapping 
of subsurface impurities. 

In the case of a tunnel contact, the oscillatory part of the conductance can be expressed 

by 

^^oc^sin2^^7, ^F7»l; (47) 
Gc 7 

where kp — ^/kp + 2meV/ h? is the wave vector of electrons that are passing through the 
orifice and Zq is the depth of the defect under the surface; 7 is the distance between the 
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contact and the defect. Comparing this to the results found for a baUistic contact, where 



oc -^27^ y^kpzo + (ppj (48) 



Gc a? {zq) 

(here is the phase the electron acquires after reflection by the contact), we see that 
although both oscillations have similar arguments, the expression for the baUistic case has 
an extra phase ip^ which depends nonlinearly on the wave vector A;, making the signal hard 
to identify. Secondly, the adiabatic condition, being an essential assumption in the ballistic 
model, cannot be readily achieved experimentally. 

Therefore, choosing the tunnel contact for experimental application seems most sensible. 
In that situation wc can expect the information in the conductance signal about a defect's 
whereabouts to be twofold: the amplitude will decrease with growing distance 7, whereas the 
frequency of the oscillation is expected to increase upon enlarging the distance from contact 
to defect. The actual experiment would consist of sensitively measuring ^{V) curves on a 
tight grid of p coordinates. The lateral positions of defects could then be identified as the 
centers of radially symmetric patterns in this signal. Next, the depth of an impurity should 
be derived from the period of the oscillation in the ^{V) curve at p^. 

Assuming the numerical parameter 2T:m*gkp/li^ = 0.5 introduced in Sec. Ill (which can 
be shown to be applicable for hard wall scatterers with atomic radius) and choosing the 
orifice to be located exactly above the defect (7 = zq), the amplitude of the oscillation is 
expected to be lO^^Gc for Zo=3 nm (with kp = 10^° m~^). 

Note that the choosing value of interaction constant is rather large. We use the such 
value of the parameter to show more clear the investigated effects in illustrations. For real 
value of parameter g ^ 10~^^erg ■ cm^, which can be estimated from an electron effective 
scattering cross section ~ lA^, the relative amplitude of oscillations is 10^^ ^ 

Comparing this to previous STS experiments [19], where signal-to-noise ratios of 5 • 10~^ 
(at 1 nA, 400 Hz sample frequency) have been achieved, we should be able to measure defects 
located more than 10 atomic layers under the surface. 

As the period of the oscillation becomes longer for small zq, the minimum discernable 
depth will be determined by the maximum voltage that can be applied over the junction. 
For example, 30 mV is sufficient for probing a quarter of a conductance oscillation caused 
by a defect at 1 nm depth. The increase of the noise level inherent to measuring at elevated 
voltages will not pose a problem, as the amplitude of the signal is much higher for small 
depths. 

Finally, the anisotropy of the electronic structure will have to be taken into account. 
Materials with an almost spherical Fermi surface such as Al or Au, realizing the condition 
of a free electron gas, are expected to be most suitable. Furthermore, deviations of spherical 
symmetry might be used as a secondary proof for the effectiveness of the method, i.e. in the 
case of Au(lll), where the 'necks' in the Fermi surface should cause a defect to be invisible 
when probed exactly from above. 

This research was supported partly by the program "Nanosystems, nanomaterials and 
nanotechnology" of National Academy of Sciences of Ukraine. Ye. S. Avotina wishes to 
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